Introduction: Aq Qala industrial state is one of the most important infrastructures in Golestan province. The layout of this area in the vicinity of Gorganrood and Gharesoo rivers connected by Gaz channel has increased its vulnerability against medium and large floods. March 2018 flood in Golestan province has just revealed many weak points of this area versus any upcoming floods. Therefore, in this paper, the author tries to remedy Aq Qala industrial state resilience against possible future floods. It is also attempted to find a reliable solution to protect this area. Various numerical models have been developed to simulate flood inundation to delineate the floodplain zones bordering the rivers and calculate the associated risk considering various return periods. Numerical models to simulate floods are categorized into (a) one-dimensional (1D) models, (b) two dimensional (2D) models, and (c) one-dimensional river flow models coupled with two-dimensional floodplain flow (1D-2D) models. In this paper, the two-dimensional behavior of flow flood around the state MIKE-HD model is used to simulate flood flow in the region. To achieve this, initially, the MIKE-HD model was calibrated and validated for the rivers and flood plain at the study area. Subsequently, the flood inundation is simulated using MIKE-HD for various return periods. The simulated flood inundation is validated using March 2018 flood. This flood was one of the biggest flood in the previous 500-year by which many houses were evacuated for three weeks, intensifying by snowmelt for another two weeks. This flood engaged several cities and villages in Golestan Province in Iran, more than 200 km2 area at Gorganrood river neighborhood. Methodology: Many factors and parameters must be considered to simulate flood inundation around the industrial state. The first requirement is topographical data which are essential to find the flood flow direction. In this study, three datasets are used, while the first one goes back to data gathered by photogrammetric operation. More than 6799 hectares are surveyed by drone in this data set with an accuracy of about 20m to 50m. The second one goes back to the ground survey around the industrial state. The last one is related to the Gorganrood topographical data provided by Golestan province Regional Water Organization. After gathering topographical data sets together, an event is considered to calibrate the flood model in the MIKE-HD software with two parameters: water levels and flood inundation. To achieve this, the flood that occurred in March 2018 is chosen as the event. In this flood, many parameters such as the observed flood plain by aerial images and water levels are considered calibration parameters. In the model, the bridge laid on Aq Qala-Gogran road over Gharehsoo river was taken into account to make the model more accurate. Finally, the model is developed and evaluated versus various floods with different return periods. Results show that the maximum water level for a 500-year return period is-15. 94m CD. Also, results show that flood flow behaviour in the area has a two-dimensional behaviour with low speed. This indicates that protection alternatives like river dikes haven’ t got any problem to be collapsed by scouring. Also, it is concluded that the obtained water levels for different return periods converge to each other in the highest value of return periods. Results and Discussion: Obtained results show that the flood flow has a two-dimensional behavior in the study area for all return periods, so floods leading themselves by water level changing and inundation. This confirms two options for industrial state protection against possible future floods. The first one is the construction of walls to protect the area while increasing the foundation levels for the structures can be the second alternative. Changing foundation levels is too complicated now because of building construction and any movement. Besides, to protect the industrial area from runoff inundation construction, some pump stations can be recommended, especially flood times. The pump station is a remedy to convey produced runoff by the site at rainy times on the back of the protection wall, especially in inundation time. Also, results show that a level of about-15. 94 CD can be a reasonable level to construct the protection wall as the ultimate level. Not that wave and wind-setup freeboard must be added to this level because of inundation enduring. Conclusion: In this study, flood flow around Aq Qala industrial state was studied. Results show that the best option for protecting the study area is a flood wall and ten pump stations at the lower level areas. Highlights about the study can be categorized as follows: 1. Simulation of March 2018 flood at the study area 2. Finding final options for the protection of the industrial state against possible future floods. 3. Flood levels for various return periods and their ultimate values 4. Flood behaviour for the flat area like Aq Qala industrial state is showing itself as level fluctuation rather than a high-velocity fluid flow 5. To remedy flood inundation at the study area, considering wave and wind set up as a freeboard to design the protection wall is essential.